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ABSTRACT 

We present an investigation of the clustering of the faint {i'ab < 24.5) field galaxy population in the redshift range 
0.2 < z < 1.2. Using 100,000 precise photometric redshifts extracted from galaxies in the four ultra-deep fields of 
the Canada-France Legacy Survey, we construct a set of volume-limited galaxy samples. We use these catalogues to 
study in detail the dependence of the amplitude A^n and slope S of the galaxy correlation function w on absolute 
Mb rest- frame luminosity, redshift, and best-fitting spectral type (or, equivalently, rest-frame colour). Our derived 
comoving correlation lengths for magnitude-limited samples are in excellent agreement with measurements made in 
spectroscopic surveys. Our principal conclusions are as follows: 1. The comoving correlation length for all galaxies 
with —19 < Mb — 51ogh < —22 declines steadily from z ~ 0.3 to z ~ 1. 2. At all redshifts and luminosity ranges, 
galaxies with redder rest-frame colours have clustering amplitudes between two and three times higher than bluer 
ones. 3. For both the red and blue galaxy populations, the clustering amplitude is invariant with redshift for bright 
galaxies ( — 19 < Mb — 5\ogh < —22). 4. At z ~ 0.5 for less luminous galaxies with Mb — 5 log ft ~ —19 we find higher 
clustering amplitudes of ~ 6h~^ Mpc. 5. The relative bias between redder and bluer rest-frame populations increases 
gradually towards fainter magnitudes. Among the principal implications of these results is that although the full galaxy 
population traces the underlying dark matter distribution quite well (and is therefore quite weakly biased), redder, 
older galaxies have clustering lengths which are almost invariant with redshift, and by z ~ 1 are quite strongly biased. 

Key words, observations: galaxies - cosmology:large-scale structure of Universe- Astronomical data bases: Surveys 



^ 1. Introduction information concerning the galaxy formation process and 

tracing the evolution of bias as a function of scale and mass 
In the cold dark matter model structures grow hierarchi- of tj^e hosting dark matter halo is one of the key objectives 
J5 . cally under the influence of gravity. Galaxies form inside of observational cosmology On large scales, (> lOh-^ Mpc) 
■ - ' "haloes" of dark matter (jWhite fc Rees, 1978). Because structure growth is largely driven by gravitation (where we 
these haloes can only form at the densest regions of the dark measure the correlations between separate haloes of dark 
matter distribution, the distribution of galaxies and dark matter); however on smaller scales (< 1 h-^) non-linear 
matter is not the same; the more stroiigly clustered galaxies effects generally associated with galaxy formation domi- 
are said to be "biased" (|KaiseH, ll98J; [Bardeen et al.|, |1986D ^ate the structure formation process. In this paper we must 
with respect to the dark matter distribution. This relation- bg^r in mind that although we measure a clustering sig- 
ship between dark and luminous matter provides important ^al to around 0.1 deg at z - 1 this corresponds to around 
~ 3h~^ Mpc and ~ 2h~^ Mpc at z ~ 0.5, which means 



Send offprint requests to: H. J. McCracken that our observations are mostly in non-linear to strongly 

* Based on observations obtained with non-linear regimes where environmental effects play an im- 

MegaPrime/MegaCam, a joint project of CFHT and portant role in the evolution of structure. 
CEA/DAPNIA, at the Canada-France-Hawaii Telescope Qn linear scal es, as theory and si mulations have 

i^^^V "^^'^^ j' operated by the National Research Council ^^^^^ (^^^ example [je^Ah^Tindl (ITolsl ) or lWeinberg et Zl 

(NRC) or Canada, the Institut JNational des Science de lUmvers /lonn,iK\ ii i j. • i-j. j r~j i j.j. j 

^r,, ^^ , ,, '. 1 J 1 r> 1, u o • 4--C rr<T^Tnc\ c (I2004I I I . the clustermg amplitude oi dark matter decreases 

or the Centre National de la Recherche Scientifaque (CNRS) oi ' fv .,, i i -n rr i • r i ; i ;i i i 

France, and the University of Hawaii. This work is based in part steadily with redshift. If galaxies perfectly traced the dark 

on data products produced at TERAPIX and the Canadian matter component, then their clustering amplitudes would 

Astronomy Data Centre as part of the Canada-France-Hawaii decrease at each redshift slice, in step with the underly- 

Telescope Legacy Survey, a collaborative project of NRC and ing dark matter. However, as the galaxy distribution is bi- 

CNRS. ased, stellar evolution intervenes to complicate this picture; 
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in effect, the actual measured clustering amplitudes are a 
complicated interplay between the underlying dark mat- 
ter component and how well the luminous matter traces 
this galaxy distribution, or how efficiently galaxies form. 
Understanding fully the evolution of galaxy clustering re- 
quires, therefore, some insights into the galaxy formation 
process. Thanks to large spectroscopic redshift surveys we 
now have a much more complete picture of the evolution of 
the g alaxy luminosity function with redshift (jllbert et al.l . 
l2005f ) and how the fracti on of galaxy ty pes evolves w i th red - 
shift (jZuccaet al.U2006D . For example. Ilibert et al.l (|2005D . 
using first-epoch data from the VVDS redshift survey have 
shown that the luminosity function brightens considerably 
between z = 0.3 and z = 1, with M* increasing by one 
or two magnitudes at z ~ 1. We must take this into ac- 
count when comparing clustering amplitudes measured at 
the same absolute luminosities in different redshift ranges. 

In the local Universe, million-galaxy redshift surveys 
have greatly expanded our knowledge of galaxy clustering 
at low redshift. We now have a broad idea how the dis- 
tribution of galax ies depends on their intrinsic luminosity 
and s pectral type (jNorberg et al.l . [20011 . l200HZehavi et all . 
|2005() . In general, these works have shown to a high preci- 
sion that at the current epoch more luminous galaxies are 
more clustered than faint ones, and that similarly redder 
objects have higher clustering amplitudes than bluer ones. 
Other works have shown th at slope of the galaxy cor relation 
depends on spectral type (jMadgwick et al.L l2003f l. These 
studies have indicated that, in general, more luminous, red- 
der, objects are more strongly clustered than bluer, fainter 
galaxies. Some studies have also related physical galaxy 
properties, such a total in ass in stars, with the clustering 
properties (jLi et al.l . l2006t) . But how do these relationships 
change with look back time? 

At intermediate redshifts (z ~ 1), however, our 
knowledge is still incomplete. Multi-object spectrographs 
mounted on ten-metre class telescopes have made it possi- 
ble to construct samples of a few thousand galaxies. The 
first studies investigating galaxy clustering as a function 
of the object's rest frame luminosity and colour for large 
galaxy samples at z ~ 1 have now appeared IjMeneux et al.l . 
120061 iPollo et all . l2006t ICoil et al.l . l2006h . Unfortunately, 
these surveys typically contain ~ 10'^ galaxies, which are 
enough to select objects either by type and absolute lumi- 
nosity, but not, for instance, to apply both cuts simultane- 
ously. 

These works confirm some of the broad trends seen 
at lower redsh i ft and with magnitude-limited samples 
(|Le Fevre et al.L l2005a| ) but are still not quite large enough 
to investigate in detail how galaxy clustering depends si- 
multaneously on more than one galaxy property. For exam- 
ple, one may investigate the dependen ce of clustering ampl i- 
tude within a volume- limited sample ([Meneux et al.l . l2006l ). 
but one may not, as yet, investigate simultaneously samples 
selected by type and absolute luminosity. Unfortunately, 
even with efficient wide-field multi-object spectrographs, 
gathering redshift samples of thousands of galaxies at red- 
shift of one or so requires a significant investment of tele- 
scope time. 

Photometric redshifts offer an exit from this impasse, 
and represent a middle ground between simple studies us- 
ing imaging data with magnitude or colour-selected sam- 
ples and spectroscopic surveys. Several attempts have been 
made in the past to carry out galaxy clustering studies 



with photometri c redshifts, mo s tly using the Hubble deep 
field data sets (^Arnouts et all l2002t iTeplitz et all . 120011 : 
Magliocchetti & Maddox, 199!|). However, such works ei- 
ther suffered from sampling and cosmic variance issues or 
poorly-controlled systematic errors. The advent of wide- 
field mosaic cameras like MegaCam (jBoulade et al.l . |2000() 
has made it feasible for the first time to construct samples 
of tens to hundreds of thousands of galaxies from z ~ 0.2 all 
the way to z '--^ 1 and beyond. Two key advances have made 
this possible; firstly, rigorous quality control of photometric 
data, and secondly, the availability of much larger, reliable 
training samples reaching to faint [j! ^ 24) magnitudes 

In this paper we will describe measurements of galaxy 
clustering derived from a large sample of galaxies with ac- 
curate photometric redshifts in Canada-France legacy sur- 
vey (CFHTLS) deep fields. These fields have been observed 
repeatedly since the start of survey operat ions in June 
2003 as part of the on-going SNLS project (jAstier et al.l 
l2006l ) and consequently each filter has very long integra- 
tion times (for r and i bands the total integration time in 
certain fields is over 100 hours). A full description of our 



photo metric redshift catalogue can be found in lllbert et al.l 
(|2006D . Containing almost 100,000 galaxies to i' < 24.5 we 
are able to divide our sample by redshift, absolute lumi- 
nosity and spectral type. These photometric redshifts have 
been calibrated using 8,000 spectra from the VIRMOS-VLT 
deep survey fVVDS: JLeFevre et al.l (|2005b[ )V In addition, 
our sample has sufficient volume to provide reliable mea- 
surements of galaxy clustering amplitudes at redshifts as 
low as z ~ 0.2; and we are thus able to follow the evo- 
lution of galaxy correlation lengths over a wide redshift 
interval. In the lower redshift bins, the extremely deep 
CFHTLS photometry means it is possible to measure clus- 
tering properties of a complete sample of objects as faint as 
Mb — 5 \og{h) ~ — 18 (at z ^ 0.2 we have large numbers of 
very faint objects with Mb — 51og(/i) ~ — 15 although we 
do not consider them here). Moreover, by using all four in- 
dependent deep fields of the Canada-France legacy survey 
we are able to robustly estimate the amplitude of cosmic 
variance for each of our samples. 

Our objective in this paper is to determine, first of all, 
how the observed properties of galaxies determines their 
clustering. We are able to carry out such an investigation 
of galaxy clustering strength for samples of galaxies selected 
independently in absolute luminosity, rest-frame colour and 
redshift. 

Our paper is organised as follows: in Section |2| we de- 
scribe how our catalogues were prepared and how we com- 
puted our photometric redshifts; in Section |3| we describe 
how we measure galaxy clustering in our data; our results 
are presented in Section |4l Finally, our discussions and con- 
clusions are presented in Section |5l In this work we divide 
the CFHTLS galaxy samples in three principal ways: first 
of all, we consider simple magnitude limited samples, di- 
vided by bins of redshift (described in Section |4T|) : next, 
at two fixed redshift ranges, we consider galaxy samples se- 
lected by absolute luminosity and type (Section 14. 2|) ; and 
lastly, at a range of redshift bins and for the same slice in 
absolute luminosity, we consider galaxies selected by type 
(presented in Section |44|) . 

Throughout the paper, we use a flat lambda cos- 
mology (r2„j = 0.3, Oa = 0.7) and we define 
h = i/o/100 km s~^ Mpc~^. Magnitudes are given in 
the AB system unless otherwise noted. 
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2. Catalogue preparation and photometric redshift 
computation 

We now describe the preparation of the photometric cata- 
logues used to derive our photometric redshifts. Although 
our input catalogue has already been released to the com- 
munity as part of the CFHTLS-T0003 release (hereafter 
"T03" ) , no extensive description of the catalogue process- 
ing has yet appeared in the literature; for completeness we 
provide a brief outline of the principal processing steps in 
this Section. 

These photometric catalogues were released by the 
TERAPIX data centre to the Canadian and French com- 
munities as part of the T03 release and have been made 
public world-wide one year later. They comprise observa- 
tions taken with the MegaCam wide-field mosaic camera 
(jBoulade et al.l . |2000[ ) at the Canada-France-Hawaii tele- 
scope between June 1st, 2003 and September 12th, 2005. 
Full details of these observations, data reductions, cat- 
alogue preparation and quality assessment steps can be 
found on the TERAPIX web page£|: however, we now out- 
line the principal steps in data reductions and catalogue 
preparation. 



2.1. Production of stacked images 

MegaCam is a wide-field CCD mosaic camera consisting of 
36 thinned EEV detectors mounted at the prime focus of 
the 3.6m Canada France Hawaii Telescope on Mauna Kea, 
Hawaii. The detectors are arranged in two banks. The nom- 
inal pixel scale at the centre of the detector is 0.186"/pixel; 
the size of each detector pixel is 13.5/i. All observations for 
the CFHTLS are taken in queue-scheduling mode. Each of 
the four fields presented in this paper have been observed in 
all five MegaCam broad-band filters primarily for the super- 
novae legacy survey. After pre-processing (bias-subtraction 
and fiat-fielding) at the CFHT, images are transferred to 
the Canadian astronomy data centre (CADC) for archiv- 
ing, and thence to TERAPIX at the lAP in Paris for pro- 
cessing. At TERAPIX, the data quality assessment tool 
"QualityFITS" is run on each image, which provides a 're- 
port card' in the form of a HTML page containing informa- 
tion on galaxy counts, stellar counts, and the point-spread 
function for each individual image. Catalogues and weight- 
maps are also generated. At this point each image is also 
visually inspected and classified. 

In the classification process galaxies are divided into 
four grades according to seeing and associated image fea- 
tures (for instance, if the telescope lost tracking or other 
artifacts were present). Only the two highest-quality grades 
are kept for subsequent analysis. 

After all images have been inspected, and bad images 
rejected, an astrometric and photometric solution is com- 
puted using the TERAPIX tool scam p which co mputes a 
solution simultaneously for all filters (Bertinl . [20061) . Finally, 
this astrometric solution is used to re-project and co-add 
all images (and weights) to produce final stacked image. 
All of these steps are managed from an web-based pipeline 
environment. The internal r.m.s. astrometric accuracy over 
the entire MegaCam field of view is always less than one 
MegaCAM pixel (0.186") 



|http : //terapix . iap . f r/rubrique . php?id_rubrique=208| 



2.2. Quality assessment 

Galaxy number counts, stellar colour-colour plots and in- 
completeness measurements have been calculated for all 
four deep stacks in all five bands. By examining the po- 
sition of the stellar locus in each field in the u — g vs 
g — r and g — r vs r — i colour-colour planes we see 
that the photometric zero-point accuracy field to-field is 
~ 0.03 or better. Detailed comparisons between CFHTLS- 
wide survey fields and overlapping Sloan Digital Sky Survey 
fields show systematic errors of a comparable amplitude. 
This degree of photometric precision is essential if we 
are to compute accurate photometric redshifts. A full 
list of the characteristics of release T03 can be found at 
http : //terapix . iap . f r/cplt/tab_t03ym . html 



2.3. Catalogue generation 

Once images have been resampled and median-combined 
for each field we use swarp to pro duce a "chi-squared" 
detection image (jSzalav et al.l I1999D based on the g', r' 
and i' stacks (the pixel scale on each image in all fields 
and colours is fixed to .186" /pixel. Next, sextractor 
([Bertin fc Arnouta . I1996D is executed in "dual-image" 
mode on all stacks using the chi-squared image as the 
detection image. This method "automatically" produces 
matched catalogues between each stack as in all cases the 
detection image remains the same. We note that, given the 
strict criterion on the image seeing used to select input 
images in the CFHTLS stacks, all deep stacks are approx- 
imately seeing-matched, with the median seeing on each 
final stack in each band of ^ 1" . This means one can safely 
use dual-mode detection. We use S extractor's map:_auto 
Kron-like "total" galaxy magnitudes (jKronl . flDSOT ) . At faint 
magnitudes, where the error on the Kron radius can be 
large, our total magnitudes revert to simple 3" diam- 
eter aperture magnitudes. After the extraction of cata- 
logues redundant information is removed from each band 
and a "fiag" column is added to the catalogues con- 
taining information about the object compactness using 
the "local" measuremen t of the object's half-light radius 
([McCracken et al.l . l2003[ ). A mask file, generated automati- 
cally and fine-tuned by hand, is used to indicate areas near 
bright stars or with lower cosmetic quality, and this infor- 
mation is incorporated in the object flag. Objects used in 
the subsequent scientific analysis are those which do not 
lie in these masked regions, are not saturated, and are not 
stars. 



2.4. Photometric redshift computation and accuracy 

A full description of our m ethod used t o com pute photo- 
metric redshifts is given in lllbert et al.l (|2006f ). Briefly we 
use a two-step optimisation process based on firstly the 
bright sample (to set the zero-points) and the full sam- 
ple (which optimises the templates). This new template set 
is then used to compute photometric redshifts in all four 
fields. In this paper we consider photometric redshifts com- 
puted using only the five CFHTLS filters {u*griz). This 
is true even in fields where additional photometric infor- 
mation is available (for example, CFHTLS-Dl field where 
there is supplementary BVRIJK photometry). This ap- 
proach was taken to ensure that field-to-field variation in 
photometric redshift accuracy as a function of redshift was 
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kept to a minimuin. Our photometri c redshifts a re essen- 
tially identical to those presented in [Ilbert et al.' with the 
exception that in the D2 field we use additional ultra-deep 
u* imaging kindly supplied by the COSMOS consortium; 
this serves to equalise the m* integration time between the 
fields. We separate stellar sources from galaxies by using a 
combination of sextractor f luxjradius parameter and the 
best fitting spectral template. 

We emphasise that a key aspect of our photometric red- 
shifts is that extensive comparisons have been made with 
large database of spectroscopic redshifts (|Le Fevre et all . 
l2005a[). In particu lar, we draw attention to Figures 9 and 
10 of lllbert et al.l which show photometric redshift accu- 
racy and the fraction of catastrophic errors as a function 
of redshift. For galaxies with i' < 24 in the redshift range 
0.2 < z < 1.2 the photometric redshift accuracy in the Dl 
field, expressed as cfazH^ + z), is always less than 0.06; 
in the redshift range 0.2 < z < 0.6 it is less than 0.04. 
Catastrophic errors are defined as the number of galaxies 
with \zs — Zp\/{1 + Zs) > 0.15 where Zg is the spectroscopic 
re dshift and Zp the photometric redshifts. From Figure 10 
in lllbert et aLl we can see that the fraction 77 of objects with 
catastrophic redshift errors is better than 5% in the redshift 
range 0.2 < 2: < 1.2 for objects with 22.5 < iab < 24.0. 

Although there are smaller numbers of spectroscopic 
redshifts in the other fields, some useful comparisons can be 
made; using 364 publ icly-availabl e spectra from the DEEPl 
project. Figure 14 in lllbert et al.l shows that the dispersion 
5z/{l + Zs) is 0.03 in the redshift interval 0.2 < z < 1.2. 
In the D2 field we have carried out an additional compar- 
ison with spectroscopic redshifts obtained by J. P. Kneib 
and collaborators in the context of the COSMOS project. 
This test, making use of 335 i' < 24 spectroscopic redshifts, 
shows that, once again, in the interval 0.2 < z < 1.2, our 
photometric redshift errors 5z/{l + Zg) are ~ 0.035. 

During the preparation of this article, an independent 
comparison has been carried out by members of the DEEP2 
team between their large spectroscopic sample and the 
CFIITLS-T03 photometric redshifts presented here. They 
find an excellent agreement between, comparable to the val- 
ues presented here for the other fields, for more than 20,000 
galaxies in the D3 survey field. 

We would like to use photometric redshifts for objects 
fainter than the Iab < 24.0 VVDS spectroscopic limit. We 
can define another figure of merit, the percentage of objects 
with o-^p(68%) < 0.15(1 + zp), where CTsp(68%) is the 68% 
phot ometric redsh ift error bar. This is plotted in Figure 
15 inl llbert et al.l and gives an indication of how good the 
photometric redshifts are beyond the spectroscopic limit. In 
the interval 0.2 < Zp < 1.5, this is always better than 80% 
for all four CFHTLS deep fields even as faint as i' < 24.5. 

2.5. Computing absolute magnitudes and types 

We measure the absolute magnitude of each galaxy in 
UBVRI standard bands {U Bessel, B and V Johnson, R 
and / Cousins). Using the photometric redshift, the asso- 
ciated best-fit template and the observed apparent mag- 
nitude in one given band, we can directly measure the k- 
correction and the absolute magnitude in any rest-frame 
band. Since at high redshifts the fc-correction depends 
strongly on the galaxy spectral energy distribution it is 
the main source of systematic error in determining abso- 
lute magnitudes. To minimise /c-correction uncertainties. 
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Fig. 1. Rest frame AIjj — Mv colour as a function of B-band 
absolute magnitude for the Dl field. Each panel from top 
to bottom corresponds to the redshift bins used in this pa- 
per. The points show the four different best-fitting spectral 
types. In the colour version of this figure (available electron- 
ically) red, magnenta, green and blue points correspond to 
Coleman et al. Ell, Sbc, Scd, and Irr templates. The right- 
hand panels show the colour distribution for each redshift 
slice. 



we derive the rest-frame luminosity at A using the object's 
apparent magnitude closer to A x (1 -I- z). We use either the 
r', i' or z' observed apparent magnitudes according to the 
redshift of the galaxy. The procedure is described in the 
Appendix A of Ilbert et al. (2005) where it shown that this 
method greatly reduces the dependency of the fc-corrections 
on galaxy templates. 

Galaxies have been classified using multi-colour infor- 
mation in a similar fashion to other works in the lit- 
erature (for example iLin et all . Il999l ; IWolf et al.l . 120031 : 
IZucca et al.l I2006D . For each galaxy the res t-frame colours 
were matched with four templates from IColeman et al.l 
(198(|) (hereafter referred to "Coleman, Wu and Weedman" 
or "CWW" templates). These four templates have been 
optimised using the V VDS spectroscopic redshifts, as de- 
scribed in lllbert et all , and are presented in Fig. 2 of this 
work. Galaxies have been divided in four types, correspond- 
ing to the optimised E/SO template (type one), early spi- 
ral template (type two), late spiral template (type three) 
and irregular template (type four). Type four includes also 
starburst galaxies. Note that, in order to avoid introducing 
dependencies on any particular model of galaxy evolution, 
we did not apply templates corrections aimed at accounting 
for colour evolution as a function of redshift. 

We show in Figure [T] the rest-frame colour distribution 
of the galaxies for each type. Type one galaxies comprise 
most of the galaxies of the red peak of the bimodal colour 
distribution. The other types are distributed in the blue 
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peak. Galaxies become smoothly bluer from type one to 
type four respectively. 



3. Measuring galaxy clustering 

3.1. Introduction 

There are two principal approaches which may be used to 
measure the clustering of objects with photometric red- 
shifts. One is simply to isolate galaxies in a certain red- 
shift range using photometric redshifts, and then to com- 
pute the projected correlation function ■w{9) for galaxies 
in this slice, as has long been done for magnitude-limited 
samples. However, the additional information provided by 
photometric redshifts on the bulk properties of our slice (its 
redshift distrib ution) allows us to use the Limber's equation 
(jLimbed . I1953D to invert the projected correlation function 
and recover spatial correlations at the effective redshift of 
the slice. These computations are easy to perform and are 
relatively insensitive to systematic errors in the photomet- 
ric redshifts as one just integrates over all galaxies in a given 
redshift slice; it has already been used extensively in smaller 
surveys and is usually the method of choice when only 
small numbers of galaxies or poorer-quality photometric 
redshifts are available, and has bee n used extens i vely o ver 
the past few y ears (see, for example iDaddi et al.l (|200lD or 
lArnouts et al.l ((199^). It has the disadvantage that it pro- 
vides only limited information on the shape of the angular 
correlation function as one measures a correlation function 
integrated over a given redshift slice. 

A second approach is to decompose the redshift of each 
galaxy into it's components perpendicular (r^) and parallel 
(tt) to the observer's line of sight, and then to compute a full 
two-dimensional correlation function ^(r^, tt) based on pair 
counts of galaxies in both directions. Finally, one computes 
the sum of this clustering amplitude in the direction paral- 
lel to the line of sight, Wp. In spectroscopic surveys, this has 
the advantage of removing the effect of redshift-space dis- 
tortions caused by infall onto bound structures. This tech- 
nique has been u sed successfully f or many spectroscopic 
redsh i ft surveys (iDavis &: Feeblest Il983t iLe Fevre et al.l . 
119961 : IZehavi et al.l . I2005D and some attempts have been 
made to apply i t to samples with lower-accuracy photo- 
metric redshifts (jPhleps et al.ll200a) . It has the advantage 
that it can provide direct information on the shape of the 
correlation function but this comes at the price of much 
greater sensitivity to systematic errors in the photometric 
redshifts (for example, integration over a much larger range 
in redshift space is necessary). We plan an analysis using 
this technique in a forthcoming article, but in this paper 
we adopt a conservative approach, as we are primarily in- 
terested in the overall clustering properties of our galaxy 
samples. 



From these masked catalogues of object positions, we 
measure w(9), the projected angular correlation function, 
using the standard Landy fc Szalavi (|1993i) estimator: 



w{9) 



DD - 2DR + RR 
RR 



(1) 



where DD, DR and RR are the number of data-data, 
data-random and random-random pairs with separations 
between 9 and 9 + 69. These pair counts are appropri- 
ately normalised; we typically generate random catalogues 
with ten times higher numbers of random points than input 
galaxies. 

An important point to consider is that, of course, the 
precision of our photometric redshifts are limited. In a given 
redshift interval, zi < z < Z2 it is certainly possible that a 
given galaxy may be in fact outside this range. To account 
for this, we employ a we i ghted estimator of w{9), as sug- 
gested bv lArnouts et al] ( 2002|) . In this scheme we weight 
each galaxy by the fraction of the galaxy's probability dis- 
tribution function enclosed by the interval zi < z < Z2. 

In this case for each pair we must now compute 



DD^^P'P';DR^ Y^ 



P' 



(2) 



where P*'-' is the integral of object's probability distri- 
bution function in the redshift interval zi < z < Z2. The 
total 'effective' number of galaxies becomes 



N,, 



E 



P' 



(3) 



We may then fit w{9) = A^9 
amplitudes and slopes A^ and S 
the "integral constraint" which arises because the mean 
density is estimated by the survey itself. 



^ to find the best-fitting 
after first correcting for 



c = ^ f fuj{9)dnidn2 



(4) 



where il is the total area subtended by our survey. We 
compute C by numerical integration of Equation (j4]), dis- 
counting pairs closer than 1" , corresponding to the resolu- 
tion limit of our data. The integral constraint correction is 
subtracted from the power law. We always fit our data on 
scales where theintegral constraint correction is negligible. 

We compute w in a series of logarithmically spaced bins 
from log{9) = -3 to log(6') = -0.2 with 6log{9) = 0.2, 
where 9 is in degrees. In the following sections we will con- 
sider both measurements derived from each field separately 
(using that field's redshift distribution) and also measure- 
ments constructed from the sum of pairs over all fields (in 
which case we use a combined, weighted redshift distribu- 
tion derived from all fields). 



3.2. Projected angular clustering 

We first use our photometric redshift catalogue to produce 
a galaxy sample generated using a given selection criterion, 
for example either by absolute or apparent magnitude, red- 
shift or type. This same selection criterion is applied to 
catalogues for all four fields. 



3.3. Derivation of spatial quantities 

We can associate each value of A^ and S with a correspond- 
ing comoving correlation length, tq by making use of the rel- 
ativistic Limber equation (Peebles, 19 80). For further dis- 
cu ssion of this me t hod se e, for example [Paddi et all (|200lh 
or lArnouts et al.l (|1999r ). If we assume that the spatial 
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correlation function can be expressed as ^(r) = {'>'/fa) '^ 
(where 7 = 1 — 5 and T is the incomplete Gamma function) 



r(7/2) 



[J{dN/dz)dz 



where dN/dz represents the redshift distribution and 9 
is the angular separation on the sky. Here g(z) depends only 
on cosmology and is given by: 



g{z) = {dx/dz)-^x^-"'F{x) 

with the metric defining x and F{x): 
ds^ = c^dt^ - a^[dxlF{xf + x2(sin2 M^^)] 

Defining 



rliz^s) 



Qg{z){dNldzfrQ{z)'<dz 
' Qg{z){dNldzYdz 



(6) 



(7) 



(8) 



and assuming that the correlation amplitude is constant 
over our redshift interval zi < z < Z2 and that w(0) = 
■Aw9^~^ we can write: 



An. 



?'o(^eff)V7r' 



r((7 - l)/2) Q g{z){dN/dzfdz 
r(7/2) [^^l{dN/dz)dz]^ 



(9) 



Thus, given a measurement A^ of the correlation func- 
tion amplitude for the redshift slice under consideration 
and a knowledge of that slice's redshift distribution we can 
derive rQ{zcs)- 

Given that we have four separate fields, we may de- 
rive To from a "global" w{9) derived from the sum of pairs 
over all fields and computed using the average, weighted 
dnj dz. Alternatively, we can compute rg for each field using 
that field's redshift distribution and individual correlation 
function amplitudes; the final value of tq is then calculated 
simply as the average over all four fields. We find that, in 
general, these methods agree for most samples. However, 
in some cases the error bars are larger with the field-to- 
field measurements. This is discussed in more detail in the 
following section. 

3.4. Error estimates on w{9) 

We have investigated two different approaches to estimate 
the errors on our measurements of rg. As described above, 
for each of our four fields, we compute rg from a fit to a 
measurement of w{9) and the redshift distribution for that 
field. We then compute the variance in rg and 7 over all 
fields. 

In the second method, for each sample, we compute the 
sum of the pairs over all fields at each angular bin. The 
error bar at each angular bin is then calculated from the 
variance in w{9) over all fields. If Wav is the mean corre- 
lation function then Wi is the correlation function for each 
field, then the error over the n fields of the CFHTLS is 
given as 



a^ = 



n 



(n-l) 



w.r 



(10) 
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Fig. 2. The amplitude of the angular correlation w as a 
function of angular separation 9 (in degrees) for 17.5 < i' < 
24 galaxies selected in the four deep fields of the CFHTLS 
in a range of redshift slices. The error bars correspond to 
the amplitude of the field-to-field variance over all fields. 
The dashed line shows the best-fitting power law correlation 
function after the subtraction of the appropriate integral 
constraint. 



where n = 4 for the CFHTLS (note that Wav is only 
used in this computation and not in any other part of the 
analysis) . 

A "global" correlation length is determined using the 
average redshift distribution over four fields. In this case, 
the error in rg is computed from the error in the best fit- 
ting j4u,. This error is, in turn, computed from the covari- 
ance matrix derived using the Levenberg-Marquardt non- 
linear fitting rout ine, as presented in Numerical Recipes 
iPresset all ((l986[ ). 

We find that in most cases the error bars in rg estimated 
by these two different methods are consistent. However at 
lower redshifts ranges (0.2 < z < 0.6), where the num- 
bers of galaxies is smaller, field-to-field dispersion is higher 
than the global errors. Further investigations reveal this is 
due to the presence of a single field (d2) which has anoma- 
lously higher correlati on lengths. Properties of this field are 
discussed in detail in iMc Cracken et al.l ()2007f ). This is un- 
doubtedly due to the presence of very large structures at 
z ~ 0.3 and z ~ 0.7 in this field We believe that our 'global' 
correlation function provides a more robust estimate of the 
total error and we adopt this measurement for the remain- 
der of the paper. 



4. Results 

4.1. Magnitude limited samples 

Can photometric redshifts be used to make reliable mea- 
surements of galaxy clustering at z ^^ 1? In this Section 
we will construct a sample similar to those already in the 
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Fig. 3. The comoving correlation length, tq as a function 
of redshift for the four combined CFHTLS fields (filled 
squares) compared to literature values (open symbols) com- 
puted for a galaxy sample limited at i' < 24.0. For these 
fits, both 7 and Vq are free parameters. The error bars on 
these fits are computed from the field-to-field variance in 
each bin in w, as described in Section [31 



literature in order to address this question. We consider 
galaxies selected by redshift and apparent magnitude, fn 
each field we divide our catalogues into a series of bins 
of width Sz = 0.2 over the range 0.2 < z < 1.2. In 
each bin, galaxies with 17.5 < iAB < 24.0 are selected 
to match the criterion used b y the VIMOS-V L T deep sur- 
vey (VVDS) as presented in iLe Fevre et al.l (|2005af) (we 
note that the Megacam instrumental i' magnitudes are 
ve ry close to the C FH12K instrumental / magnitudes used 
in iLe Fevre et al.[ ). Following the procedures outlined in 
Section [31 we measure the weighted pair counts at each 
angular separation for each field and sum them together. 
Equation (|T]) is used to derive a "global" w{d). The error 
bar at each bin of angular separation is computed from the 
variance of the individual measurements of ^1(6) over all 
four fields, as described in Equation (fTUj) . These results are 
illustrated in Figure [2l which shows the amplitude of w for 
all four fields for a range of redshift slices. We note that in 
all redshift slices except the lowest-redshift one, at inter- 
mediate scales, w{9) is well represented by a power law. In 
addition our error bars are reassuringly small. This global 
w(0) is then fitted with the usual power law, correcting 
for an integral constraint corresponding to a total area of 
3.2deg2. 

In calculating the correlation amplitude ro at the ef- 
fective redshift of each slice we use a redshift distribution 
derived from the weighted, summed dn/dz from each field. 
Our results are displayed in Figures [31 and [31 They are con- 
sistent with th e measurements from th e VVDS deep survey 
(open svmbols: lLe Fevre et al.l . l2005al) . which was based on 
a much smaller sample of ~ 7000 spectroscopic redshifts. 



Fig. 4. The best-fitting slope 7, as a function of sample me- 
dian redshift for the four CFHTLS fields combined (filled 
squares) compared to literature measurements (open sym- 
bols) for a sample selected with 17.5 < i' < 24.0 at each 
redshift slice. 



Our data does show some evidence for a decline in the cor- 
relation amplitude strength in the interval 0.5 < z < 1.1, 
as well as a slightly higher slope, in contrast with this ear- 
lier work. This effect will be discussed in more detail in 
Section [5l 

4.2. Luminosity dependent clustering at z ^ 0.5 and z ^ 1.0 

The main difficulty in interpreting Figures [31 and [4l is that 
at each redshift slice the sample's median absolute lumi- 
nosity changes significantly as a consequence of the selec- 
tion in apparent magnitude. This is evident if one considers 
Figure [5l which shows a two-dimensional image of the ob- 
jects distribution in the absolute magnitude-redshift plane. 
Slices extracted at lower redshift are dominated by galaxies 
of intrinsically low absolute luminosity. 

To investigate the dependence of galaxy clustering am- 
plitude on absolute luminosity, in this Section we extract 
samples in two fixed redsliift intervals selected from the 
absolute-magnitude/luminosity plane. We consider galax- 
ies between 0.2 < z < 0.6 and 0.7 < z < 1.1. In each 
redshift range we select a minimum absolute magnitude 
(shown by the solid lines in Figure [5]) so that the median 
redshift of galaxies selected in each slice of absolute lumi- 
nosity is approximately constant. At each redshift interval 
we separate the galaxy population into "early" (types one 
and two) and "late" (types three and four). We also con- 
sider samples comprising all galaxy types. These samples 
are illustrated in Figure [6l where the median rest-frame 
Johnson [B — I) colour is plotted as a function of median 
rest-frame Mb magnitude. In both high- and low- redshift 
slices, changes in the absolute magnitude bin produces the 
largest changes in rest-frame colours. It is also clear that 
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Fig. 5. Gray-scaled histogram showing the distribution of 
galaxies as a function of absolute magnitude and redshift for 
four CFHTLS fields for all galaxy types and for an apparent 
magnitude limit of i' < 24.5. 
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Fig. 6. Median rest-frame (B-I) colour versus median rest 
frame S— band absolute luminosity for samples at 0.2 < 
z < 0.6 (open symbols, sohd hnes) and 0.7 < z < 1.1 
(filled symbols, dotted lines). For each redshift range we 
show all galaxy types (circles) , types one and two (squares) 
and types three and four (triangles). 



galaxy populations become progressively bluer at higher 
redshifts (dotted lines and filled symbols), at the same bin 
in absolute magnitude. 

Our type-selected correlation functions for galaxies in 
the redshift range 0.2 < z < 0.6 are displayed in Figure [7l 
and for the higher redshift range in Figure [S] These plots 
show the amplitude of the angular correlation function w 
as a function of angular separation, 6, for different slices 
of absolute magnitude. In each panel we show correlation 
functions measured for the red and blue (early and late) 



populations. The size of the error bars at each angular sep- 
aration corresponds to the amplitude of the field-to-field 
cosmic variance computed over the four fields. The long- 
dashed and dashed lines show the fitted amplitudes for 
the red and blue populations respectively. For the higher 
redshift bin, we superimpose the fitted amplitudes at the 
same bin in absolute luminosity at the lower redshift inter- 
val. At all redshifts and absolute luminosity ranges, galaxies 
with redder rest-frame colours are more clustered than their 
bluer counterparts. 

Because of the strong covariance between 7 and tq it is 
useful to consider contours of constant x^ a-t each slice in 
absolute magnitude. This is shown in FigurelHl The vertical 
and horizontal lines mark an arbitrary reference point. Solid 
lines indicate galaxies with 0.2 < z < 0.6 and dotted lines 
those with 0.7 < z < 1.1. From these plots we see a gradual 
increase in comoving correlation length as a function of 
absolute rest-frame luminosity. We also see some evidence 
for a slight decrease in the comoving correlation length at 
given fixed absolute luminosity between higher and lower 
redshift slices. 

Figures [10] and [TT] show the comoving correlation func- 
tion x^ contours for red and blue populations in both red- 
shift ranges; and correspond to the angular correlation func- 
tions presented in Figures [7] and Figures [H Figure [9] shows 
the comoving correlation length for the full galaxy popula- 
tion in both redshift ranges. 

The results presented in this Section are summarised in 
Figures [T^l [131 and in Tables [1] and [D These figures show the 
best-fitting correlation amplitude as a function of absolute 
luminosity for the three different samples (early, late and 
the full galaxy population) in the two redshift ranges (0.2 < 
z < 0.6 and 0.7 < z < 1.1) considered in this Section. 

Considering these plots, several features are apparent. 
Firstly, at all absolute magnitude slices and in both redshift 
ranges, early- type galaxies are always more strongly clus- 
tered (higher values of ro) than late-type galaxies. Secondly, 
we note that clustering amplitude for the late-type popula- 
tion is remarkably constant, remaining fixed at r^ 2h~^ Mpc 
over a large range of absolute magnitudes and redshifts. The 
behaviour of the early- type population is more complicated. 
For the 0.2 < z < 0.6 bin, we some evidence that as the 
median luminosity increases, the clustering amplitude of 
this population decreases, from around ^ 6/1"^ Mpc for the 
faintest bins, to ^ 5h~^ Mpc. We note that the difference 
in clustering amplitude between the early and late popula- 
tions is smaller for the higher-redshift bin. We also note that 
the clustering amplitudes we derive for our blue and full- 
field galaxy populations are cons iderably lower than those 
reported by iNorberg et al.l ()2002[ ) at lower redshifts. 

We carried out several test to verify the robustness of 
the higher correlation amplitudes observed for the fainter 
red galaxy population at z ^ 0.5. Selecting galaxies with 
redder rest-frame colours (classified as type one) in larger 
bins of absolute magnitude we also measure higher cluster- 
ing amplitudes for objects fainter than AIb — 5 log h ^ —19. 
We also note that the origin of the large error bar for the bin 
at Mb — 5 log ft. = — 19 is due to the presence of structures 
in one of the four fields; interestingly, for the Dl field, for 
this faint red population, the correlation function does not 
follow a normal power-law shape. Our resulting error bars 
reflect this behaviour, but it is clear that for certain galaxy 
populations, for instance the bright elliptical population, 
simple power-law fits are not appropriate. 
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Fig. 7. The amplitude of the angular correlation w a,s a. function of angular separation 9 (in degrees) for i' < 24.5 
galaxies in the rcdshift range 0.2 < z < 0.6. Each panel shows galaxies selected in a different, independent, slice in 
absolute luminosity. Squares represent the early-type population (types 1 and 2) whereas triangles show the late-type 
population. The dashed and long-dashed lines shows the adopted best- fit. The error bars at each angular separation 
correspond to field-to-field variance measured over the four survey fields. 



In the redshift range 0.7 < 2: < 1.1 we have compared 
our measurements of pure luminosity dependent cluster- 
ing (i.e., without type selection) to works in the literature 
comput ed using snialler s ampl es of spectroscopi c redshifts, 
namelv 'PoU o et all (|2006D and lCoil et"all (|2006[ ). shown in 
Figure [T^ as the open triangles and open stars respectively. 
Their points should be compared with the full circles de- 
rived from our measurements. In general the agreement is 
acceptable, although for higher luminosity bins, our ampli- 
tudes are below the measurements from the DEEP2 survey 
(although it seems that the amplitude of their error bars is 
perhaps underestimated) . 



For all the plots previously shown in this Section we 
fitted simultaneously for the slope and amplitude of the 
galaxy correlation function. In Figure [14] and Figure [15] 
we summarise our results from our low and high redshift 
samples. Figure [15] shows the slope 1 — 7 as a function of 
absolute luminosity for early-type, late type, and full-field 
populations at high redshifts. Figure [Ml presents the results 
from the 0.2 < z < 0.6 sample. Error bars are computed 
from the field-to-field variance. 

Interestingly, we find for the higher redshift bin (0.7 < 
z < 1.1) the slope is relatively insensitive to absolute mag- 
nitude. However, at lower redshifts, luminous red galaxies 
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Fig. 8. Similar to Figure [T] the amplitude of the angular correlation w as a function of angular separation 9 (in degrees) 
for red and blue i' < 24.5 galaxies in the redshift range 0.7 < z < 1.1. The dashed and long-dashed lines shows the 
adopted best-fit for the same bin in absolute luminosity in the lower redshift range 0.2 < z < 0.6. The error bars at each 
angular separation correspond to field-to- field variance measured over the four survey fields. 



have a steeper correlation function slope than fainter galax- 
ies. A similar effect is observed in the SDSS and two- degree 
field surveys at lower redshifts (JNorberg et al.l . l2002f ) . 



amplitude of the pure type four population and find that in 
this case clustering amplitudes are lower than the combined 
sample of types three and four.) Error bars are computed 
from the field-to-field variance. 



4.3. Clustering of bright early type galaxies 

In Figure \W\ we plot the clustering amplitudes of bright 
early-type galaxies in our survey; filled squares indicate 
galaxies of types one and two, and open squares represent 
a pure type one sample. As expected, the clustering ampli- 
tudes of the pure type one population (with overall redder 
rest-frame colours) are higher than the combined type one 
and two samples. (We have also measured the clustering 



Several authors have presented clustering measurements 
as a function of cither absolute lum inosity, type or redshift. 
For example. iMeneux et al.l (|2006f ) described measurements 
in the VVDS spectroscopic redshift survey of the projected 
correlation function for early- and late-type galaxies. Their 
galaxies are classified in the same way as in this paper, 
using CWW templates. However, in their sample galaxies 
were selected by apparent magnitude; at z ~ 1, their rest 
frame luminosities are comparable to the brightest galaxies 
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Fig. 9. The comoving correlation length rg and slope 6 — 
1 — 7 as a function of median absolute luminosity and red- 
shift for low-redshift (0.2 < z < 0.6; solid lines) and high 
redshift (0.7 < z < 1.1; dotted hues) populations for four 
fields of the CFHTLS-deep. Each panel shows contours of 
constant chi-squared values for 1,2, and 3a confidence lev- 
els with the plus symbol marking the minimum chi-squared 
value. For convenience, lines of constant slope and comov- 
ing correlation length are plotted. 



Fig. 11. The comoving correlation length tq and slope 
(5 = 1 — 7 as a function of absolute luminosity and redshift 
for late- type populations for the four fields of the CFHTLS. 
Each panel shows contours of constant chi-squared values 
for 1,2, and 3a confidence levels with the plus symbol mark- 
ing the minimum. The solid and dotted lines shows galaxies 
in the redshift bin 0.2 < z < 0.6 and 0.7 < z < 1.1 respec- 
tively for the same range of absolute luminosities. 
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Fig. 10. The comoving correlation length tq and slope 
6 = 1 — 7 as a functions of absolute luminosity and redshift 
for early- type galaxies for the four fields of the CFHTLS. 
Each panel shows contours of constant chi-squared values 
for 1,2, and 3a confidence levels with the plus symbol mark- 
ing the minimum. The solid and dotted lines shows galax- 
ies in the redshift bin 0.2 < z < 0.6 and 0.7 < z < 1.1 
respectively for the same range of absolute luminosities. 
We note that in the -18.5 < Mb - 5log{h) < -17.5 and 
— 19.5 < Mb — 5log{h) < —18.5 luminosity ranges the clus- 
tering amplitude of early type galaxies at different redshifts 
are well separated. 



in our sample. We compare these z ^ 1 with our data; they 
are shown as the open circles in Figure 1161 Finally, the 
open triangles represent measurements from iBrown et al.l 
(|2003t ) who measured clustering of red galaxies selected us- 
ing three-band photometric redshifts in the NOAO wide 



survey. Their results are above ours by at least one or two 
standard deviations. In general we note that our results 
are lower than literature measurements and speculate that 
this could be the consequence of a slight loss of signal due 
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Fig. 12. The comoving correlation length ro as a function of 
median absolute luminosity and type for objects in the red- 
shift range 0.2 < z < 0.6. Filled circles show the full galaxy 
population. In addition to type selection, the galaxy sam- 
ple is selected in one-magnitude bins of absolute luminos- 
ity. Triangles and squares represent the late and early-type 
populations respectively. 
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Fig. 13. Similar to Figure [T^ here we show measurements 
from our 0.6 < z < 1.1 sample. Points from the literature 
from the measurements made by the DEEP 2 and VVDS 
surveys. 
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Table 1. Low redshift sample (0.2 < z < 0.6). Columns 
show the median rest-frame i?— band absolute luminosity, 
total number of galaxies over the four fields, median ab- 
solute rest-frame B-I colour, effective redshift, best fitting 
correlation length and slope. 



our use of photometric redshifts. The broad trend seen in 
our measurements is that the clustering amplitude of bright 
early- type galaxies does not change with redshift. 

4.4. Redshift-dependent clustering 

How does the comoving correlation length of the various 
galaxy populations investigated here depend on redshift? 
Over the full range redshift range (0.2 < z < 1.1), as is 



Fig. 14. Best- fitting slope 1 — 7 of w as a function of median 
absolute luminosity at 0.2 < z < 0.6 for the early- type 
population (filled circles), late- type population (triangles) 
and full population (open circles). Slopes are plotted as a 
function of the median absolute magnitude in each slice of 
one magnitude in width. 
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Table 2. High redshift sample (0.7 < z < 1.1). Columns 
show the median rest-frame _B— band absolute luminosity, 
total number of galaxies over the four fields, median ab- 
solute rest-frame B-I colour, effective redshift, best fitting 
correlation length and slope. 



apparent from Figure [SI this measurement is only possi- 
ble for the brightest galaxies; at higher redshifts, intrinsi- 
cally fainter galaxies drop out of our survey. For each of 
the redshift bins used in Figure [5] we selected galaxies with 
—22.0 < Mb — 5\ogh < —19 and measured ro and 7 as 
above. The derived amplitudes for this sample are shown 
as the stars in Figure [T71 We also selected at each red- 



McCracken et al.: Galaxy clustering in the CFHTLS 



13 





, , 1 , , , , 1 , 








i' < 24.5 




_ 




* Late-type 




- 


1.2 


— • Early— type 




— 




O All 




r 


i 1 


- 




J 


j „ : 




( : 




1 


) 






T 




( 


) 


0.8 







---. 




-■ 




I 

1 1 1 1 1 1 1 


1 


J 
1 1 1 1 1 1 


1 1 1 1 




-19 -20 


-21 


-22 






M 


B-5 


log 


(h) 





o 
ft 



2 - 



-22. 0< MB-51ogh<-l£ 
i' < 24.5 



• All types 
■ Type 1+2 



"■ i ^ 



A Type 3+4 



0.5 



Fig. 15. As in Figure [T^] but for the redshift range 0.7 < 
z < 1.1. Slopes are plotted as a function of the median 
absolute magnitude in each bin of one magnitude in width. 



Fig. 17. Redshift dependence of comoving galaxy correla- 
tion length tq for a series of volume limited samples for 
early types (squares) , late types (triangles) and for the full 
sample (filled circles). 
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Fig. 16. Clustering amplitude of luminous red galaxies. 
Open and filled squares show measurements for type one 
and type one and two combined galaxy samples. Other 
points show measurements from the literature for early type 
galaxies selected using a variety of methods. 



shift bin early type galaxies (types 1 and 2) and late type 
galaxies (types 3 and 4); these are represented by squares 
and triangles in Figure [T71 We find that the median ab- 
solute magnitude at each bin is Mb — 51og/i ~ —19.6 for 
the early types and Mb — 5 log h ^ 19.4 for the late types. 



The full galaxy population has an absolute magnitude of 
Mb - 5logh ^ -19.5. From OA < z < 1.2, these values 
changes by at most 0.1 magnitudes. As in previous plots, 
the size of the error bars represent cosmic variance errors 
over the four fields. 

We find that in the redshift range 0.4 < z < 1.2 probed 
by our survey, early-type galaxies are always more clus- 
tered than late type galaxies as we have already found in in 
the previous Sections. Moreover, the difference in clustering 
amplitudes between these two populations is approximately 
constant with redshift. (Note that the galaxy samples ex- 
amined here correspond to essentially the brightest bins 
plotted in Figures fT2l and Figure [13]). We also find that the 
clustering of the —22 < Mb — 5\ogh < — 19 luminosity- 
limited full galaxy sample (i.e., without including a type 
selection) decreases steadily from z ^ OA to z ^ 1.2. 

4.5. Relative bias 

We can also compute the relative bias between different 
galaxy populations at different redshifts. At each redshift 
range in Section [Ol fO. 2 < z < 0.6 and 0.7 < z < 1.1) we 
compute the relative bias b as follows: 

b = ba/h^as{a)/as{b) (11) 

We adopt the usual definition for as (jPeeblesl Il980t) , 



CT8 



^iC^{ro/m 



Where is C7 is a constant which depends on 7: 
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Fig. 18. Top panel: relative bias between early and late- 
type populations at z ~ 0.5 (open circles) and z ^ I (filled 
circle s). The open t r iangle s an d star show measurem ents 
from ISwanson et all (|2007t) and iMarinoni et all (|2005D re- 
spectively. Bottom panel: relative bias between the low- and 
high-redshift full galaxy sample. For both panels, biases 
are plotted as a function of the median absolute rest- frame 
magnitude in slices of one magnitude in width. 



In Figure [T51 we plot the relative bias between the early 
and late-type populations for our low and high redshift sam- 
ples (open and filled circles respectively) as a function of 
absolute rest-frame luminosity. Error bars are computed 
from the field-to-field variance. From our data it is clear 
that the relative bias between the early and late type pop- 
ulations declines between z ~ 0.5 and z ~ 1. 

We compare our relative bias between blue and red 
galaxy populat ions with measurements in the literature. 
IMarinoni et al.l (2005) divided their spectroscopic sample 
into a red one with rest-frame {B — I) > 0.95 and a 
bluer one with [B — I) < 0.68. By measuring the prob- 
ability distribution function (PDF) of these two popula- 
tions, they were able to be measure the relative bias of 
relatively bright {Mb — 51og(/i) ~ —21) galaxies in the 
interval 0.6 < z < 1.1, indicated by the starred point in 
Figure [TH] Recentlv. ISwanson et al.l (|2007( ) investigated the 
relative bias between red and blue spectroscopic galaxy 
samples in the Sloan Digital Sky Survey as a function of 
absolute rest-frame red magnitude. Their results are repre- 
sented in Figure [T51 as open triangles (with an approximate 
offset applied to convert to rest-frame blue magnitudes used 
in this paper). 

From Figure [TS] we see that o ur measurements at 
z ^ 1 agree with ISwanson et al.l and IMarinoni et al.l . 
ISwanson et al.l findings are similar to ours: they observe an 
increase in the relative bias between early and late types at 
fainter magnitudes. However, at z '^ 0.5, our relative bias 
measurements are considerably higher than either measure- 
ments z '^ or z '^ 1. 



In this paper we have used a sample of 100,000 photometric 
redshifts in the CFHTLS legacy survey deep fields to inves- 
tigate the dependency of galaxy clustering on rest-frame 
colour, luminosity and redshift. The first sample we con- 
sidered comprises a series of magnitude-selected cuts sam- 
pling the full galaxy population from 0.2 < z < 1.2. We 
find that the clustering amplitude decreases gradually from 
3/i"i Mpc at z ~ 0.3 dechning to 2h-^ at z - 1.0. The de- 
clining correlation amplitude for the full galaxy population 
at least to z ~ 1 indicates that the field galaxy population 
must be weakly biased, as this trend follows that of the 
underlying correlation amplitudes of the dark matter. 

We next repeated the same experiment (measuring 
galaxy clustering in narrow redshift bins) and imposed an 
additional selection by absolute luminosity and rest-frame 
colour. Selecting galaxies by slices of absolute luminosity 
the steady decline in comoving correlation length found for 
samples selected in apparent magnitude is even more pro- 
nounced (Figure [T7|) . The luminous field galaxy population, 
dominated by blue star-forming galaxies at z ~ 1, is clearly 
only weakly biased with respect to the dark matter distri- 
bution. 

Turning to rest-frame colour-selected samples at all red- 
shift ranges we consistently find that galaxies with red- 
der rest-frame colours are more strongly clustered than 
those with bluer rest-frame colours (Figure |17|) . Such an 
effect has long been observed for galaxie s in the local 
Univers e (for exampl e Norberg et al., '2002'; 'Ze havi et al.l 
l2005; Lovedav et al.l . 1999; Guzzo et al., 199'$~~and~at 
higher redshifts for samp les selected by ty pe and luminos- 
ity (jMeneux et all . 120061 : [Coil et all . [20061 ) . Numerical sim- 
ulatio ns find a similar effect: For example, [Weinberg et al.l 
()2004D show that older, redder galaxies are more strongly 
clustered. This is a generic prediction from most semi- 
analytic models and hydrodynamic simulations of galaxy 
formation: older, more massive galaxies formed in regions 
which collapsed early in the history of the Universe. At the 
present day such regions are biased with respect to the dark 
matter distribution. 

For the brightest ellipticals (-22 < Afs -5 log h < -19) 
in our survey, we find that their clustering amplitude does 
not change with redshift (Figure \T7\ . indicating that at 
z '--^ 1 the elliptical population must be strongly biased 
with respect to the underlying dark matter distribution. 
Comparing our measurements for objects with redder rest- 
frame colours with those of other surveys, we find similar 
clustering amplitudes. Reassuringly, as we demonstrated 
in Section 14.3! sub-samples of galaxies with redder rest- 
frame colours produce even higher correlation amplitudes 
(Figure dH). 

In a second set of selections we considered the depen- 
dence of galaxy clustering on luminosity and type in two 
broad redshift bins: 0.2 < z < 0.6 and 0.7 < z < 1.1 
(we leave a 'gap' in the range 0.6 < z < 0.7 to ensure 
that there is no contamination between high and low red- 
shift ranges). Once again, for the most luminous objects 
{Mb — 5 log h ~ —20) the correlation amplitude is approxi- 
mately constant between these two redshift bins. However, 
for fainter red objects, at a fixed absolute luminosity, we 
see a decline in correlation amplitude between z ~ 0.4 
and z ~ 1; the same is true for samples selected purely 
by absolute magnitude. We find no evidence for a change 



McCracken et al.: Galaxy clustering in the CFHTLS 



15 



in clustering amplitude at the same luminosity for the blue 
population with redshift. 

At 0.2 < z < 0.6, where we are complete to Mb — 
5 log h < — 17, we find that red galaxies with Mg — 5 log h ~ 
—20 are more strongly clustered than bluer galaxies of the 
same luminosity. Moreover as the sample rest frame lumi- 
nosity decreases to Mb — 5logh — 18 the clustering am- 
plitude rises from ~ 4/i~^ Mpc to ~ 6h^^ Mpc. A simi- 
lar effect has been re ported in larger , low-redshift samples 
in th e local universe (jSwanson et al.l . l2007t iNorberg et al.l . 
I2002D . where both the Sloan and 2dF surveys have found 
higher clustering amplitudes for redder objects fainter than 
L* . Some evidence for this effect has also b een reported in 
numerical simulations (jCroton et al.l . [2007l ). which indicate 
that this behaviour arises because faint red objects exist 
primarily as satellite galaxies in halos of massive, strongly 
clustered red galaxies. This means that less luminous, red- 
der objects reside primarily in higher density environments 
at z ~ 0.5. This is in agreement with recent studies of 
galaxy clusters at intermediate redshift which indicate a 
rapid build-up of low luminosity red galaxies in clusters 
since z ~ 1 Ivan der Wei et al.l (|2007[ ). Our survey is not 
deep enough to probe to equivalent luminosities at z ~ 1. 

Conversely for the redshift bin at 0.7 < z < 1.1 we see 
that for the full galaxy population more luminous objects 
are more strongly clustered: ~ 2h^^ Mpc for galaxies with 
Mb— 5 log h ^ —19 and ~ 4:h~^ Mpc for galaxies with Mb— 
51og/i ~ —21. At all luminosity bins, galaxies with redder 
rest-frame colours are always more strongly clustered than 
bluer galaxies. 

In both redshift ranges, we measured the slopes of the 
correlation function as a function of redshift, luminosity 
and rest-frame colour. At z '--^ 1 we observe that redder 
galaxies have steeper slopes; at lower redshifts however dif- 
ferent galaxy populations have identical slopes. At these 
redshifts, we find steeper slopes in our most luminous bin; 
at higher reds hifts, no su c h obv ious trend is apparent (in 
contrast with iPollo et all (2006') , who saw a clear depen- 
dence of slope on absolute luminosity) . 

We have also computed the relative bias between red 
and blue galaxies at z ~ 1 and Z'-^^O.S. Atz'^l our re- 
sults agree with measurements in the literature. Our mea- 
surements at z '--^ 0.5 are significantly above measurements 
made at z ~ 0. Interestingly, our results show that the 
relative bias between early and late types increases gradu- 
ally for samples selected with fainter intrinsic luminosities, 
which is consistent with the results presented for our inves- 
tigation of galaxy clustering. 

Concluding, we may summarise our results as follows: 
firstly, for samples of galaxies with similar absolute lu- 
minosities, galaxies with redder rest-frame colours are al- 
ways more strongly clustered than their bluer counterparts. 
Secondly, for the bluer galaxy populations, the correla- 
tion length depends only weakly on absolute luminosity. At 
lower redshifts, we find some evidence that redder galax- 
ies with lower absolute luminosities are more strongly clus- 
tered. For the entire galaxy population (red and blue types 
combined) we find that as the median absolute magnitude 
increases, the overall clustering amplitude increases. For 
our the most luminous red and blue objects, the clustering 
amplitude does not change with redshift. 

The overall picture we draw from these observations is 
that the clustering properties of the blue population is re- 
markably invariant with redshift and intrinsic luminosity. 



In general, galaxies with bluer rest-frame colours, which 
comprise the majority of galaxies in our survey, have lower 
clustering amplitudes (typically, ^ 2h^^ Mpc) than the 
redder populations. The clustering amplitude of the blue 
population depends only weakly on redshift and luminos- 
ity. This is consistent with a picture in which bluer galaxy 
types exist primarily in lower density environments. 

In contrast, the clustering amplitude of the low- 
luminosity red population is lower at higher redshifts. In 
FigurefTTlwe see that for the luminous (Mb — 5 log h ~ —20) 
red population, the correlation amplitude does not change 
with redshift. Moreover, at a fixed absolute luminosity, the 
correlation amplitude of the full galaxy population and 
the magnitude-selected galaxy population decreases from 
z '^ 0.4 to z ~ 1.1, in step with the underlying dark matter 
distribution. 



6. Conclusions 

We have presented an investigation of the clustering of 
the faint {i' < 24.5) field galaxy population in the red- 
shift range 0.2 < z < 1.2. Using 100,000 precise photo- 
metric redshifts extracted in the four ultra-deep fields of 
the Canada-France Legacy Survey, we construct a series of 
volume-limited galaxy samples and use these to study in 
detail the dependence of the amplitude A_w and slope S 
of the galaxy correlation function w on absolute Mb rest- 
frame luminosity, redshift, and best-fitting spectral type 
(or, equivalently, rest-frame colour). Our principal conclu- 
sions are as follows: 

1. The comoving correlation length for all galaxies with 
— 19 < Mb — 5 log ft, < — 22 declines steadily from z ~ 0.3 
to z ^ 1. 

2. At all redshifts and luminosity ranges, galaxies with 
redder rest-frame colours have clustering amplitudes be- 
tween two and three times higher than bluer ones. 

3. For both the red and blue galaxy populations, the 
clustering amplitude is invariant with redshift for bright 
galaxies with -22 < Mb - 5\ogh < -19. 

4. At z ~ 0.5, less luminous galaxies with M^- 5 log ft ~ 
— 19, we find higher clustering amplitudes of ^ 6ft^^ Mpc. 

5. The relative bias between populations of redder and 
bluer rest-frame populations increases gradually towards 
fainter magnitudes. 

The main implications of these results is that although 
the full bright galaxy population traces the underlying dark 
matter distribution quite well (and is therefore quite weakly 
biased), redder, older galaxies have clustering lengths which 
are almost invariant with redshift, and must therefore by 
z ~ 1 be quite strongly biased. In addition, at z ~ 0.5 
there is some evidence that fainter red objects are more 
strongly clustered than ^ L* galaxies at these redshifts. 
This is consistent with a picture in which fainter red objects 
exist primarily as satellite galaxies in galaxy clusters. 

It is tempting to interpret our results in terms of stud- 
ies which show that the number density of massive, lu- 
mino us galaxies evolves little from z ^ 1 to the present 
dav ^Cimatti et al.1 . l2006t IZucca et all . l2006t ICowie et all 
|1996[ ). In our survey, the clustering amplitudes of bright el- 
lipticals are already 'fixed in' at z ~ 1. Most of the changes 
in the clustering amplitude occur in the fainter galaxy pop- 
ulation. However, a full understanding of the processes at 
work here will require mass-selected samples covering a 
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larger interval in redshift. Such samples will become pos- 
sible in the near future with the addition of near-infrared 
data to the CFHTLS survey fields. 
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